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Abstract 

> 

' We consider the following singularly perturbed elliptic problem 

' e'^Aii — u + ii^ — 0, u>0 in fl, —- — on dO,, 

Ifi , an 

' where SI is a bounded domain in with smooth boundary, e is a small parameter, n denotes 

fT^ ' the inward normal of dQ and the exponent p > 1. Let F be a hypersurface intersecting dQ in 

the right angle along its boundary dV and satisfying a non-degenerate condition. We establish 

the existence of a solution concentrating along a surface F close to F, exponentially small 

in e at any positive distance from the surface F, provided e is small and away from certain 

critical numbers. The concentrating surface F will collapse to F as e — >■ 0. 
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1 Introduction 

We consider the following problem 

du 

e'^Au-u + uP = 0, u>0 mil and — = on 917, (1.1) 

an 

where is a bounded domain in with smooth boundary, e is a small parameter, n denotes the 
inward normal of dfl and the exponent p > 1. 
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1.1 Background and assumptions 

Problem (II. ip is known as a stationary equation of Keller-Segel system in chemotaxis [H]. It can 
also be viewed as a limiting stationary equation of Gierer-Meinhardt system in biological pattern 
formation 14 . Problem (|l.ip has been studied extensively in recent years. See [30] for backgrounds 
and references. 

In the pioneering papers, under the condition that p is subcritical, i.e., 1 < p < j^^^ when 
TV > 3 and I <p < +oo when = 2, Lin, Ni and TakagiflQ]. Ni and Takagi[3T]-[32]. established 
the existence of a least-energy solution Ue of (11.11) and showed that, for e sufficiently small, Ue has 
only one local maximum point E dfl. Moreover, H{P^) — > max i/(P) as e — > 0, where H{-) is 
the mean curvature of d^. Such solutions are called boundary spike-layers. 

Since then, many papers investigated further the solutions of (jl.ip concentrating at one or 
multiple points of Vl. (These solutions are called spike-layers.) A general principle is that the 
location of interior spikes is determined by the distance function from the boundary. We refer the 
reader to the articles [3], [7], [9], [15], [38], and references therein. On the other hand, boundary 
spikes are related to the mean curvature of dil. This aspect is discussed in the papers [4], [6], [8], 
[16] ■ [T5] . [37] J s-nd references therein. A good review of the subject up to 2004 can be found in 

m- 

The question of constructing higher-dimensional concentration sets has been investigated only 
in recent years. It has been conjectured in [30] that for any 1 < fc < — 1, problem (|1.1|) has 
a solution C/g which concentrates on a fc-dimensional subset of fi. We mention some results that 
support such a conjecture. 

In [3S] and [3S], Malchiodi and Montenegro proved that for N > 2, there exists a sequence of 
numbers -> such that problem (|l.ip has a solution f/^^ which concentrates at the whole bound- 
ary dfl (or any component of dil). In |23l[24], Malchiodi showed the concentration phenomena 
for (|l.ip along a closed non-degenerate geodesic(i.e. codimension 2) of dft in three-dimensional 
smooth bounded domain ft. Mahmoudi and Malchiodi in [21] proved a full general concentration 
of solutions along /c-dimensional (1 < k < N — I) non-degenerate minimal submanifolds of the 
boundary for iV > 3 and 1 < p < ^5^- 

In the above papers [25]- [55], the higher dimensional concentration set is on the boundary. A 
natural question is if there are solutions with high dimensional concentration set inside the domain. 
For two dimensional case, the authors [331 HD] consider problem (|l.ip with solutions concentrating 
on curves near a nondegenerate line F' connecting the boundary of at right angle. The meaning 
of nondegeneracy of the line F' can be defined similarly as the nondegeneracy of the hypersurface 
F in the sequel. The reader can also refer to the survey paper by J. Wei [3S] . 
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The main objective of the present paper is to extend the resuh in 39J and show the existence of 
concentration on interior hypersurfaces touching the boundary of f2 C M'^ . We make the foUowing 
assumptions and notations. The reader can also refer to the references [TTJ[51], as well as fT^ for 
some basic geometric results. 

(Al): Our candidate hypersurface F G 51 is a minimal surface that satisfies the following assump- 
tions: r is smooth and embedded in R'^, and intersects dfl in the right angle along its boundary 
9r = r n dil, which is a simple close curve in M'^. Let be the Laplace-Beltrami operator on F 
and ki, k2 the principal curvatures of F. As a submanifold of M'^, then we define the norm of the 
second fundamental form of F by 

\Ar\^ = kl + kl 

and the mean curvature of F by A: = fci + fc2(with scalar 2). Note that k = along F. We recall 
again that n is the inward unit normal vector on dfl, and hence, it locally is also the unit normal 
vector of the curve dT because of the perpendicularity between F and dfl. Note that we here 
assume that the orientation of dT is induced from that of dfl. Since a curve on the surface dfl is 
a geodesic if and only if its normal vector is parallel to the normal vector n of dil. Therefore, by 
defining 

(/Tl 

^(2^) = -(tT''^)' y^^dvcdn. (1.2) 

we know that — / is the curvature of the geodesic on dfl passing through y € dT in the normal 
direction I'd/) of F. 

(A2): Let t be the normal of the curve dT, which is also the restriction of n on dT, I is denoted in 
(|4.3p . By defining an eigenvalue problem, which will play an important role in our considerations, 

A^f + \Ar\^f = Xf inF, df/dT + If = on dT, (L3) 

we say that F is non-degenerate if the eigenvalue problem (jl.Sp does not have zero eigenvalues. 
The reader can refer to 34 for some explanations and examples. 

1.2 The profile function w 

Before stating the main result, we introduce two functions w and Z. Let w be the unique (even) 
solution of 

w" -w + wP = 0andw>0m M, w'{0) = 0, w{±oo) = 0. (1.4) 
It is well known that the associated linearized eigenvalue problem, 

h" - h+pwP-^h^ Xh inR, [h^ = l, h e H\R) (1.5) 
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possesses a unique positive eigenvalue Aq with a unique even and positive eigenfunction Z . This 
follows for instance from the analysis in [32 . In fact, we have 



w{x) = Cp I exp 



{p-l)x 



exp 



Z 



It is easy to see that for 1 



P+i 



wP+\ Ao = -(p-l)(p + 3). 



2a 



2pCp 
p-1 



-p\x\ 



0{e 



^{2p-l}\x\^ 



Z(x) -Cpe-fP+i'l^l - ^k±^e-^P\^\ + 0(e-(3p-i)|.|)^ 



(1.6) 
(1.7) 

(1.8) 
(1.9) 

(1.10) 



where 



Cri 



ip+l) 



iv + i) 



p+1 
p-1 



f +1 dx 



1.3 Main Theorem 

Our main theorem can be stated as the following: 

Theorem 1.1. Assume that the minimal hypersurface T satisfies the nondegeneracy condition 
kl.S^) . There exists a sequence of small parameters {ei}i such that problem 11. 1]) has a positive 
solution u^, still denoting £/ hy e, concentrating along a surface T near T . Near T , takes the 
form 

I dist(y, r) \ 



(1.11) 



where w denotes the unique positive solution of problem Moreover, there exists a positive 

number cq such that satisfies globally, 

Ue{y) < exp[ -Co e^^dist{y,f) ], 



and the surface T will collapse toT as £ — > 0. 



□ 



To explain in a few words the difficulties we have encountered, let us assume for the moment 
that = K. X r is an infinite strip. In terms of the stretched coordinates (s, z) = e~^(r, y) the 
equation would look near the surface approximately like 



{s,z) e © 
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The effects of curvature and of the boundary conditions are here neglected. The hnearization of 
this problem around the profile wl-s) becomes 

Let Pi, uji, i = 1,2, ■■■ , denote the eigenvalues and eigenfunctions of —A'", then functions of the 
form 

(j>l = Ws{s)uji{sz), (ff = Z{s)uji{ez), 

are eigenfunctions associated to eigenvalues respectively —piS^ and Ao — pie^ . Many of these 
numbers are small and thus "near non-invertibility" of the linear operator occurs. These two 
effects, combined in principle orthogonally because of the L^-orthogonality of Z and Ws, are actually 
coupled through the smaller order terms neglected. 

In [TJ[20l[33], related singularly perturbed problems, the translation effect 4>l have been suc- 
cessfully treated through successive improvements of the approximation and fine spectral analysis 
of the actual linearized operator. In [251 [26] resonance phenomena similar to the "0^-effect" has 
been faced in the Neumann problem involving whole boundary concentration. In j22j-|24). |10j . 
|35) . |39|-|40|. the concentration on a fc— dimensional minimal submanifold, involving both (f)} and 
(j>f effects, has been treated via arbitrary high order approximations. 

To prove Theorem II. 1) not only the same difficulties as that in the two dimensional case in 
[39)-|40j are encountered but also more obstruction appears. In other words, in the present paper 
we need to handle more delicate resonance phenomena due to the existence of 0f -effect in higher 
dimension, as well as the strong interaction between the concentration set and the boundary due 
to the homogeneous boundary condition in Some words are in order to explain the methods 

to handle these difficulties. 

As we have stated in the above, we first neglect the boundary condition. By the suitable 
rescaling of the variables (s, z) = £~^{r, y) in a type of Fermi coordinates(cf. Lemma l2.2|) . we then 
try the inner approximate solution to the problem roughly in the form, 

v{s, z) = w(s - !(ez)) + ee{ez)Z{s - f{ez)) + (p{s, z), (1.12) 

where Z is defined in (|1.5I) . / and e are left as parameters, while (f{s, z) is L^-orthogonal for each 
z both to Ws{s — f{sz)) and Z{s — f{ez)). Solving first in (p a natural projected problem, the 
standard reduction procedure will implies that the resolution of the full problem will be reduced 
to a nonlinear, nonlocal second order system of differential equations in (/, e). The nonlinear 
system is not solvable due to the (p^ effect in the case of dimension iV = 3. This shows that 
the approximation in ()1.12|) doesn't work well as N becomes large. Hence we must improve our 
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approximation by the recurrence method used in [3T], [3S]. The principle is: the better the 

approximation, higher the chances of a correct inversion of the fuU problem to obtain a contraction 
mapping formulation of the nonlinear, nonlocal second order differential equations. To do that, we 
try the following form as our new approximation, (see [26] ) 

3 

v{s, z) = w[s - f{ez)) + ee{ez)Z(s - J{ez)) + ^ e^Vi{s, z). 

1=1 

The aim of adding the term X]f=i ^VK'^i is to cancel the error term till order 0(£^) such that 
our approximation is good enough. After very tedious but necessary computations we find that 
such (pi^s must satisfy some differential equations(cf. (I3.52[) . (I3.53p . (13. 54^ ). whose solvability relies 
on the parameter /. So we need to improve our approximation further by subtle adjusting the 
location of the concentration set, namely we take the following form of inner approximate solution, 
(see [22], [35]) 

3 2 

V = w{x) + ee{ez)Z{x) + '^^e'' ipi{x,ez) with x = s — '^^e'' fi[ez). (1-13) 

1=1 1=0 

To fulfill the objective we need to analyze the corresponding linear problem very carefully and 
conduct lots of computations. This was done in Section [3j Note that the authors in [3F also 
met the same situations for the concentration phenomena on high dimensional hypersurface of 
inhomogeneous Schrodiger equation on any N dimensional space. 

Second, we know that the inner approximation in general does not satisfy the boundary con- 
dition in (|l.ip . Note that, for the two dimensional case of (jl.ip . by using the condition that 
the limit location of concentration set, say F' as before, connect the boundary dfl orthogonally, 
the authors [39]- [40] used a type of local coordinates (previously used by M. Kowalczyk in [20] ) 
in a neighborhood of V to stretch the boundary dQ and decompose the interaction between the 
concentration set and boundary dQ in such a way that they can improve the approximation to 
satisfy the boundary condition up to order of 0(£^). Unfortunately, the coordinate system can 
not be extended to the present case. Here, in a neighborhood of the hypersurface F, we will set 
up another type of local coordinates in Lemma [131 used by K. Sakamoto in [37), called modified 
Fermi Coordinates in this paper for the convenience of notion, such that we can add boundary 
corrections to the inner expansion in (|1.13p and finally get a good local approximation. The reader 
can refer to [4] for more details on the construction of boundary correction layers. 

In the rest paper we carry out the program outlined above, which leads to the complete proof 
of Theorem II. II by the infinite dimensional reduction method introduced in IIOJ. The organization 
of the paper is as follows. In Section [2] we introduce a local modified Fermi coordinates and then 
set up local formulation of problem ([1.1|) by suitable rescaling. Sections [3] and [4] are devoted to the 
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construction of a good approximate solution to solve the problem up to order 0{e'^). Indeed, we 
first ignore the boundary condition in (jl.ip and construct an inner approximation solution by the 
induction method in Section [31 In Section 4, we add the boundary correction layers to the inner 
expansion and get the final approximations involving unknown functions (/2,e)- In Section [SJ a 
gluing procedure from [1(1 reduces the nonlinear problem to a projected problem on the infinite 
strip 6, while in Section [S] and [71 we show that the projected problem has a unique solution for 
the pair of (/2,e) in a chosen region. The final step is to adjust the parameters /2 and e which 
is equivalent to solving a nonlocal, nonlinear coupled second order system of differential equations 
for the pair (/2,e) with boundary conditions. This is done in Sections [8l and [9l 

Finally, we remark that our arguments can be easily extended to show the existence of inte- 
rior concentration phenomena approaching minimal non-degenerate hypersurfaces and possessing 
interaction with boundary for problem (jl.ip on bounded domain of general dimension, see Re- 
marks [321 and [121 However, for the clearness of presentation, we here only consider the case with 
concentration on hypersurfaces in three dimension. 



2 Local Formulation of Problems 



The main objective of this subsection is to set up a local suitable coordinate system near the 
minimal hypersurface F and then express problem (jl.ll) in local form. 



2.1 Geometric notions 



Let 7o(-) iU— >Fcribea smooth parameterization of the surface F, where F is the minimal 
surface in Section 1 and D := {y e : \y\ < 1} is the unit disk. For further computational 
convenience, we will use an isothermal representation 70 : B — > F. Namely, we use 70 that satisfies 

2 





2 


970 






dy^ 



■2f . / 970 970 \ „ 



(2.1) 



Therefore, the tangent vectors d-fo/dy^ and d-fo/dy^ have the same length fi(y) > and are 
mutually orthogonal. 

Recall the notations given in subsection ll.il For our future setting-up of problem (|l.ip . here we 
provide a type of local coordinates modified from the standard Fermi coordinates, called modified 
Fermi coordinates for convenience of notions. This coordinate system was previously used by 
K. Sakamoto to describe the transition layer for Allen-Cahn equation in [Mj . 
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Lemma 2.1. There exist constants (5 > 0, tq > 0, which depend only on T and dfl, and a smooth 
diffemorphism 



7(.,.) : (~ro,ro) xB5^f]J°, 
where and fij" are given in fOj) and p77p suc/i i/ia< 

i. 7(0, y) = 7o(y) /or y e B, 7(r,y) e dn for y e 9D, \r\ < tq; 

|l(0,y) ^^.(y) /or yeD; 
3. '){r^y) has the following expansion, as r ^ 

l{r,y) = 70 (y) + ri^iy) + —qi{y) + -^92 (y) + 0{r^ 
where qi{y) and (72 (y) «re vector functions orthogonal to v{y); 



(2.2) 



ye 



(2.3) 



^. // iwe write 7 as 7(r, ??, p) in terms of the coordinate (r, p) m i2.13\) . then the derivative 
along the inward unit normal vector n of dfl is expressed as 



— = ^(9''- +9''— +9''-), 
dn ^(^33 \ dr dd dpP 

where at (r, 19,0), we have the expression 
g'Hr,d)^-r 



(2.4) 



.23 



.33 



(r, ??) = 2r 



(r,i?) 



970 


-2 


(91, 


d^o \ 






dp 1 


970 


-2 


970 


-\ 














d-fo - 



0(r' 



(2.5) 



dp 



2r 



dp \ dp' dp/ 



Proof. The proof can be found in [31]. For the convenience of readers, we also provide the details 
here. We extend 70 smoothly to 



(2.6) 



for some fixed constant (5 > 0. The extension is still denoted by 70 and its image by Ts- Let 
v{y) e M'^ be the unit normal of Ts at 70 (y) G T^. We now define a neighborhood VQ' of by the 
Fermi coordinates 



{y e : y = 7o(tj) + ry{y), \r\ <ro,ye Bs}. 
Now we choose ro > in (|2.7p small enough so that 

(riS" nr!)n{ 70 (y) + r^y) : |r|<ro, |y| = 1 + ^ } = 



(2.7) 
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When we deal with the portion of Sfi in f2^°, we use the coordinate £ dB x [0,(5) C D to 

parameterize a smaU region nearby dV in F, where the points (i?, 0) is sent to the boundary of 
surface dT by 70 and there holds 

We define the Fermi coordinates in (|2.7p by 7 : (— ro,ro) x Dg 1 — ^ R-^, i.e. 

i{r,y) ^io{y) + rv{y)- 

Note that this coordinate system is not suitable for the boundary of Q,. Whence the main objective 
of this proof is to make a suitable modification of 7. 
Denote by S the preimage of VQ' n dVt: 

S = T^{9rs° 

Since dO. _L F by (Al), we have 

5±aD({0}xD). (2.9) 
We also denote by C the preimage of 17^" n il: 

and by C(r) the r— slice of C: 

C(r) -{j/eDa : (r,2/) eC} for |r| < tq. 

For later use, we set 

Cs^T\^7). Cs{T)^{yens: {r,y)(,Cs} for |r| < tq, 

Since dVl and 70 are smooth, C(r) is a smooth domain, diffeomorphic to C(0) = D. Therefore, 
there exists a family of smooth diffeomorphisms 

y(r, •) : Bi — >C{r) 

parametrized smoothly by r <E (— tq, tq). Thanks to (|2.9p . we can choose Y so that 

dY 

no,2/) = y, ^(0,y) = 0, 2/eD. 
Furthermore, we make a smooth extension 

y(r,-) : ^^C,-(r) 
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such that 

dY 

Y{Q,y) = y, (0,y) = 0, y G B^, 

or (2.10) 

y(r,y) = Y{r,y), {r,y) G (-ro,ro) x D. 
Let us now define the desired diffeomorphism 7 by, called modified Fermi coordinates 

l{r, y) = l{r, Y{r, y)) = io{Y{r, y)) + rv{Y{r, y)) (2.11) 

for (r, y) G (— ?'o,''o) x E'a'- It is now straightforward to verify that 7 in (|2.1ip satisfies Lemma 
I2.1f l). By elementary computations and (|2.10l) . we find that 

^{Q,y) = v{y)^v{Y{Q,y)), 



l^iy) = Q^i^^y) = E 9FJ^(0'y) ^ ^(y)' (2.12) 
92(2/) = ^(0,2/) = >] (^^15I!r(0'2/) + 3^^^(0,2/)) ^ 



d^-f,^ , ^ 970 a^yj 
J=l 

which will show the validity of the statements (2) and (3). 

To prove (4) of Lemma [^TTl we use the coordinates (i?,p) introduced in the above. Recall that 
(i?, p = 0) parametrizes and p is chosen so that 

For y G D near i9D, we express 7(r, y) by 

7(r,y) =7(r,z?,p). (2.13) 

We also denote by n(r, 1)) the unit inward normal vector of dfl at 7(r, i), 0). Note that at p = 0(i.e. 
on dfl n fi^" ) , vectors 

dj 57 dj^-^s 
dr^ dd^ dp ' 

constitute a basis for . Hence n(r, i?) is expressed as 

where c > 0. Since jf^ , spans the tangent space of at 7(7", i?, 0), we have 

(|j>n)-0, (|j,n)=0, (n,n)=l at p = 0. (2.15) 

From (|2.14p and (|2.15p . we easily obtain 

^13 ^23 ^33 



^ ^ 3 3 ^^/^33 ^^^/^33 



10 



where the formulas depend on the inverse of metric matrix at p = 0, i.e. {g^-') = idij) ^- More 
precisely, the metric coefficients gij^s are given by, 



522 
933 



/97 97\ 
\dp' dpi 



Therefore, -S- is given by 



d_ _ 1 / 
9n , Irfi'i V 



;13 



.912 = .921 
ffl3 = 331 
523 = 532 



/57 57\ 
\(9i?' dpi' 



^33 



^ dr ^ d-d ^ dp 



Let us now take Taylor expansions of g-'''{r, -d, 0) in r at r = 0. From the expansion of 7(r, i?, 0) 
in Lemma l2. 1^ 3). we have, at p = 



(2.16) 



^7 _ Sto 
dp dp 



,2 

2" 



(i?, 0) + r^^id, 0) + 0) + 0(r3). 



By using the orthogonalities ^ J- , J- , qi J- and ^ J- together with 1^1^ . we 



find that, at p = 



^ 511 512 513 



dy 



dp 



521 522 523 
y 531 532 533 J 



If^l^y 



( {q,,^} (gi,^) \ 



(2.17) 



ai9 

970 \ 



V (91' if) 



-2M 



-2A^ 



0(? 



where 



-L 
-2M 



/djQ du\ 
\ d-& ' d-d/' 
/djo dvx /djo dv\ 
\ d-&' dp/ \di9' dp/ 



\ dp ' dp / ' 
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Therefore, we obtain 



' .9" 




gl3 














~g32 


g33 





I, 







I 870 I 
I 9i? I 





1^1 ^('fll 
I d-d I 






I dp I 



^70 



I '^70 I 



(91 



970 \ 
dp I 



I 870 I 870 \ 

1 1 1 870 ' ~^ 



I ^70 I ^/„, 97o\ \ 

I dp I ap / 



970 I ^ I 970 I 2 
9p 



9 P7o 070 

2 



^70 I **jY 



+ 0(r2). 

This completes the proof of Lemma 12.11 
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We can use the modified Fermi coordinate system in (|2.1ip to express the Laplacian operator 
A in (|l.ip locally in terms of (r, y) G (— ro,ro) x B. Since y = 'j{r,y), the standard metric in 
ri^" C is pulled back to g^kir^y): 



9jk 



0,1,2, 



(2.18) 



V dy^ dy^ / 

where if^ stands for r. Hence, the Laplace- Beltrami operator and gradient operator are defined in 
local coordinates by 

1 



A = 



1 ' 



(2.19) 



Vdetg 

where h is any smooth function and the coefficients g'^ are the entries of the inverse matrix of 
.9 — (ffij)- There also hold similar expressions for the Laplace- Beltrami operator and gradient 
operator on F. 

For later use, we also recall the Weyl's asymptotic formula, referring for example to [5 , or to 
[T7] and [2^ for further details. Let pi, LUi,i = 1,2, ■ ■ ■ , denote the eigenvalues and eigenfunctions 
of —A^ (ordered to be non-decreasing in i and counted with the multiplicity), then we have that 



Pi ^ — T/T^\ — as I — 00, 



(2.20) 



Vol(F) 

where Vol(F) is the volume of (F,^) and g is induced from 5, C is a constant depending only on 
the dimension — 1 . 



2.2 Local formulation of the scaled problem 

If we set u{z) = ■u{ez), then problem ()1.1[) is equivalent to the scaled problem 



AgU — u + u'' — in rip 



du 



= on dfl^ 



(2.21) 
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After rescaling, we denote il^ — ^/s and also = T/e. 

To get local form of problem (|2.2ip and construct the approximation to a solution of (|2.2ip . 
which concentrates near Fg, after rescaling, we also introduce the scaled modified Fermi coordinates 
in the neighborhood of F^ by 

(s,z) = (s,zi,Z2)e (-^, x©„ (2.22) 

where Be := {z e : \z\ < l/e} for some fixed constant e > 0. More precisely, we denote 

-fe{s,z) = -f{es,ez)/e. (2.23) 

However, for the portion of dfle in 0.^, we use the coordinate {0,1]) £ x [0,(5/e) C to 
parameterize a small region of dT^ in F^, where the points (9, 0) is sent to the boundary of surface 
iSFg by 7o and there holds 



As a direct consequence of Lemma 12. 1[ there holds 
Lemma 2.2. There exist a constant tq > 0, which depend only on F and dCl, such that 

1. 7,(0, z) = 7o(e^)/e for z e B„ 7,(5, z) G dn, for z e 9B„ |s| < ^; 

2. ^iO,z) = iyiez)forzeB,; 

3. 7^(5, z) has the following expansion, as s ^ 

ES^ ^ 

le{s, z) = jo{£z)/e + sv{ez) + —qi{ez) + ——q2{ez) + 0(£^s'*), z G B^, 

2 D 



(2.24) 



(2.25) 



where qi{£z) and q2{sz) are vector functions defined in i2.12\) . which are orthogonal to v{£z). 

4- If we write 7^ as ^^{s,9,vi) in terms of the coordinate {s,6,ri), then the derivative along the 
inward unit normal vector of dfte is expressed as 



dn^ y/g^ ^ ds do dij 
where at (s,0,O), we have the expression 

9'\s,6 

9'\s,6 

9'Hs,6 



(2.26) 



-es 





-2 




dp 








-2 


970 






dp 



2es 
970 -2 



+ Oie's') 



(2.27) 



dp 



2es 



djo 



dp \ dp ^ dp/ 
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Proof. The reader can refer to Lemma [2. II for the details of the proof. □ 
Obviously, in the scaled modified Fermi coordinates, it is of importance to express Laplace- 
Beltrami operator in (j2.21l) in local form. In fact, we obtain the form by careful calculations 



A= ^ + A^^ + Bo{s,ez) + Bi{s,ez) + B2{s,ez), (2.28) 



where 



Bois, ez) = ek— - e'^s(kf{ez) + kl{ez)) — + e^b{ez)s'^— , 
OS OS OS 

B,{s,ez) = -Vl^^-2sVl;^[-^), 



2 d ^ 9^ 

B2{s,ez) = e^^ai{es,ez)- \- e^^h{es,ez)——— + B^is^ez) . 

i—1 i—1 

In the above, we have denoted 

i=i ^ 
The functions b[ez)^ ai{es, ez),bi{es, ez),i = 1,2 satisfy: 

b{ez), ai{es,ez),bi{es,ez) < C(l + |es|^). 

Note that 0(e)-term in Bo{s,ez) is actually absent because F is minimal (cf. (Al)). Here we 
have denoted A'"^ the Laplace-Beltrami operator on F^. The vector qi is defined in p.25p and the 
differential operator i?3(s,ez) is of size O(e^). Whence, by using of Lemma [2.21 we can get local 
form of problem (I2.2ip . This will be done in more details in Sections |3] and ID 



3 Inner Approximate Solutions 

In this section, we neglect the boundary condition in (|2.2ip and then find a local inner approximate 
solution to solve the first equation in (|2.2ip up to order of 0{e'^). In fact, by the scaled modified 
Fermi coordinates in Lemma 12. 2[ we write down the local form of the first equation in problem 
()2.2ip and then extend it to a differential equation with unknown functions defined on the infinite 
strip & in R'^ with notations 

6 = { : X e R, z e FJ, de = { [x, z) : x e R, z e dT^}. (3.29) 

3.1 Inner formulation of the scaled problem 

As we have mentioned in subsection 12.21 in terms of the local coordinate system (s, z)(cf. (|2.22|) ). 
in the neighborhood of F^, the differential equation in (|2.2ip is locally expressed as 

Uss + A^-u - u + uP + B{u) = 0, (3.30) 
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where the Hnear differential operator B{u) is defined by 

B{u) ^ Bo{u) + Bi{u) + B^iu), 

by the notations 

Bo{u) = -e^ s{kl{sz) + kl{ez))^ + e^b{ez)s'^^, 

OS as 

B,{u)^-^^^^Ju)-2s^l^J^), 

B2{u) = e'^^a,{es,ez)^ + s'^^bi{es,ez)-^-^ + B^i^u). 

i—l * i—1 * 

We assume that, in the (s, z) coordinates, the location of concentration of the solution is 
characterized by the surface 

2 

f, :s = ^e7.M. (3.31) 

1=0 

Remark 3.1. The smooth functions /q, /i are to be determined in Sections\^ and^ More pre- 
cisely, with the help of the validity of the nondegeneracy of T in subsection \l.l[ we will choose /o 
by solving the equation fiS. 59\) with boundary condition and also the equation \3. 63]) with 

boundary condition ^■17\ l for fi . In fact, the nondegeneracy of T implies that /o is identically 
zero. While the unknown parameter f2 is to be chosen by a type of reduction procedure, which is 
equivalent to solving a system of differential equations in Section\^ (cf. i9.1]) - l97^ ). □ 

In the sequel, we always assume that /2 satisfies the uniform constraint 

II/2IU = ll/2||L-(r) + l|V^/2||L-(r) + ||Ar/2||L,(r) < ei (3.32) 

We consider a further changing of variables and define a new function v{x,z) as follows 

2 

u{s, z) = v{x, z) with X = s — '^^e^ fi{ez), z — z. (3.33) 

1=0 

We now want to express the problem in the new coordinates. Whence we need the following 
formulas 

2 



i=0 1=0 i=0 

2 2 



j=0 i=0 

dzi dyiJ dzi' dz^ds dyj dzidx' 
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Locally, this gives that u solves p.30p if and only if the function v defined in (|3.33p solves the 
following problem 

S{v) = v^^ + A^-v -v + yP + B{v) = 0. (3.34) 

Note that we will consider problem p.34|) on the whole ©(cf. ()3.29|) '). In the above we have denoted 
the linear operator 

B{v)=Bi{v)+B5{v) + BQ{v). (3.35) 
The linear operator ,64 (w), i?5(u) andi?6(w) can be expressed explicitly by 

1=0 1=0 i=0 

2 2 ^ 

4=0 i=0 

2 2 2 



1=0 1=0 i=0 



i=l j=0 ' i=l j=0 ' 



3.2 Inner approximate solutions 

In the subsection, we want to use coordinates (x, z) defined in p.33p to construct a suitable 
approximation to a solution expressed in the form, 

3 

V = w{x) + ee{£z)Z{x) + ^ s'^'Piix, ez), (3.36) 

i=l 

where w and Z are two functions given by (11.41) and (|1.5p . In the above expression, we have denoted 
(^i, i = 1,2,3, smooth bounded functions to be determined in the sequel. As we have mentioned, 
the unknown parameters /2(cf. p.3ip ) and e will be chosen in the last section by solving a system 
of differential equations(cf. (|9.ip - (|9.4p ). In all what follows, we shall assume the validity of the 
following uniform constraints on the parameter e 

l|e||6 = ||e|lL~(r) + £|| V^e||L,(r) + e'|| A^e|lL,(r) < . (3.37) 

For simplicity of notations, define 

r = I (/2,e)| the functions /2 and e satisfy p.32p and p.37p respectively |. (3.38) 

Now the key point is to choose suitable correction terms tpi, (p2, fs, and then prove that the 
approximate solution V solve problem (I3.34p up to order 0{e'^). 
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Formally, we have 



I w 2 \w J \w J \ w J 



W I 



(3.39) 



Setting = e(eZ + + ^ eVi ^^^d separating the powers of e, we get 

'i=2 



U;il+---+j3 



«=1 '=0 Jl,...j3,Ey.=' 

Whence, using elementary calculation, we collect the powers of e up to order 4 in the last formula, 
and then get the estimate 

3 4 

(w + eeZ + E '^Vi)^ - E ^' E 



■Ho 



=0 Ji,...j3,IZ iji=' 







e|5i 






w \ J 



More precisely, using p.39p . we make a decomposition 

'E-' E 



'=0 ii,---i3,Ii;iii=' 



i=3 



= ^1 + ^2 + ^3j 

where we have denoted 



1 4 



i=3 



In the above, we have also denoted that 



;=3 



where for every I — 3, 4, the component is independent of the terms (^3. 
Putting V into p.34p and expanding formally, we derive that 

3 



B{w) + e{e^A^e + Xoe)Z + 

3 3 

+ B{eeZ) + E e'A^'^, + E e'B{^,) + + = 0. 



(3.40) 



(3.41) 



(3.42) 



(3.43) 



(3.44) 



i=l 



1=1 
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3.3 Inner errors 

In this subsection, we accept that ipi^s as known functions for the moment. Then we compute all 
error terms in (j3.44p and make decomposition of all components into suitable forms according to 
the order of e. 

First of all, we calculate the term 

B{w)^B4{w)+B5iw)+Be{w). (3.45) 

Direct calculation gives that 

4 4 6 

2 — 1 i—1 2—5 

with expressions defined by 
Si^Si = 0, 

S2 = -A^foWx - (fci + kl)foWx, 

52 = {'^^ hTwxx - {kj + kl)xwx + ^fo'^^JaWxx, 

53 = -A^fiwx - {kl + kDhwx + &(/o + x^)wx, 

S3 = '^fo'^^JlWxx + + 2bfoXWx, 

Si = -A^f2Wx - (fci + kl)f2Wx, 

Si - {V^fl)^Wxx + 2fo'^^gj2Wxx,+2hV^JlWxx + 2f2V^JoWxx + 2bfiXWx. 

Note that for any i — 1, ■ ■ ■ , 4, 5'i is an odd function in the variable a;, while Si is an even function 
in the variable x. On the other hand, for any i = 5,6, the high order term Si is combination of 
powers of the parameters of /o, /i, /2 and their derivatives with smooth bounded coefficients. 
At the meantime, the linear operator Bq{w) comes from Bq{v) and can be expressed explicitly 

by 

4 4 

i=3 1=3 

where the terms B5(/o,/i,/2) and Bi(/o,--- ,/i-3), Bi(/o,--- ,/;-3),« = 3,4 are combination 
of powers of the parameters /o,/i,/2 and their derivatives with smooth bounded coefficients. 
Moreover, Bi(/o, • • • , fi-3), « = 3, 4 are odd functions in the variable x, while Bi(/o, • • • , fi-3), i = 

3. 4 are even functions in the variable x. As a conclusion, we get 

4 4 6 4 4 

B{w) = e'S, + e'S, + ^ e'S, + ^ + ^ e*B, + e^Bg. (3.46) 

2—1 2—1 2—5 2—3 2 — 3 
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Second, we compute the error 

sie'^A^eZ + XoeZ) + B{eeZ), (3.47) 
with B{eeZ) = Bi{eeZ) + B^{eeZ) + BeieeZ). There also holds 

3 3 6 

B^ieeZ) + Br^ieeZ) + Be{eeZ) = ^ e% + ^ e% + ^ e'T,. 

i=l i=l i=4 

In the above, we have denoted the following forms 
Ti = Ti = 0, T2 = T2 = 0, 

n = -A^/oeZ, - foikl + kl)eZ, - V^/oV^eZ, - 2/oV^^eZ„ 

^3 = (V^/o)2eZ,, - [kl + kl)exZ, - V^^eZ + - 2V^^ea;Z,, 

r,; = 6i(/o,/i,/2,e), 1 = 4,5,6. 

In the above, for any i = 4, 5, 6, the term 6i(/o, /i, /2, e) is combination of powers of the parameter 
/o,/i,/2,e and their derivatives with smooth bounded coefficients. Moreover, for any i — 1,2,3, 
Ti is an odd function in the variable x, while T,; is an even function in the variable x. 
In summary, we have that 

4 4 6 4 4 

s{v) = ^'^^ + E + E + E + E 

i—1 i—1 i— 5 i—3 i—3 

3 3 6 



i=l i=l i=i 

3 

■ ISF eZ + eXoeZ + Y^'\'P^,xx - +pwP-^ip,] 

i=l 

3 3 

Y e'A^^ V. + E ^'-^(^^) + ^2 + ^3. 



(3.48) 



i=l i=l 

Now, we shall write the error terms involving correction terms (pi, ip^ in a suitable form. In 
the next subsection, for any given i = 1, ...,3, we will choose ipi as the form 

aii{ez)bii{x) + ai2{ez)bi2{x), 

for some generic smooth functions a^i, 0^2, bn (odd) and 6^2 (even). Moreover, the terms an and 0^2 

do not depend on the unknown parameters /o, /i and/2. The reader can refer to p.60p and p.64p . 

Whence, we make a decomposition as 

3 44 
^£*Ai^-(^, = ^eWj + ^eWj + eWs. (3.49) 

i=l j=3 i=3 
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where 

Uj = A/'j(¥'i, • ■ ■ , Vj-2ja, • • ■ , /i-3), 
A^s =7V5((pi,--- ,f3,fo,--- ,/fc-2). 

Moreover, A/} is an odd function in the variable a;, while Afj is an even functions in the variable a;. 
From the definition of the operator B in p. 341) , we also write 

3 4 4 10 

1=1 j=3 j=3 j=5 

where, for any j = 3, 4, the components Gj and Gj do not depend on the correction terms 
<P3 and the unknown parameters /j_2, /2- In other words, we have 

Gj = Gjiy^i, ••■,</'i-2,/o, •••,/j-3), 
Gj = Gj{ipi, ...,ipj^2, fo, ■■■Jj-s), 

Gj = Gj{ipi,...,(f3,fo,---,f2)- 

Moreover, Gj is an odd function in the variable a;, while Gj is an even function in the variable x. 

For later use, using p. 421) and p.43l) . we decompose Hs into even parts and odd parts, and 
then write ^2 + V-s as 

% + H3 = ^e^pip - l)wP-^eZ + 

4 4 (3.50) 

+ pip - l)wP-^{eZ + £V»~i + E + 

1=3 j=3 

For j — 3, 4, the components J)j and S)j are independent of the terms ^j-i, <P3, 
i.e. 

2)j =2)j(^i,...,(^j_2), =I)j((^i,...,^j_2), (3.51) 

Moreover, J)^ is an odd function in the variable a;, while S)j is an even function in the variable x. 

3.4 Determinations of the inner correction terms (/Pj's 

In this subsection, by a recurrence procedure, we will choose suitable parameters /o, fi so that we 
can really find the correction terms ip2, and then improve the approximation. In fact, it will be 
shown that ipi is identically zero. 

This can be done in the following way. It is worth mentioning that the term 

e^A^eZ + eXoeZ 
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lies in the approximate kernel of the linearized problem of problem p.34p at V. We ignore this 
term for the moment and then cancel other components of the error in p.48|) with order of e lower 
than 4 by choosing suitable correction terms ipi, ■■■,(p3. Whence, for given ?/ e F, we then consider 
the problems 



'Pi,xx - + Pw'' Vi = in R, 
(^i(±oo) = 0, (fiWxdx = Q; 



(3.52) 



(^2(±00)==0, / (p2Wx'dx = Q 



(3.53) 



and 



(3.54) 



fPd^xx - + pw^ ^(P3^ -S3- S3-p{p-l)wP {eZ + ipi)ip2 - C3 inM, 

ip3{±oo) =0, / (P3W3: dx = 0. 
Jr 

In the above, we have denoted C3 — C3 +C3, with the odd part and even part given by 

C3 = B3 + Ts + AAa + G3 + 2)3, C3 = B3 + Ts + A/'a + G3 + 1)3. (3.55) 

Using ()3.52p . it is easy to show 1^1 = and we finish the first step. To proceed the second step 
and cancel the error terms of order O(e^) for the improvement of the approximation, we should 
choose the correction term (p2 by solving problem p.53p . For this purpose, first, we collect all 
terms of order 0{e'^) in 5(V), which has the form £^^2 with 

^2 = ^2 + ^2 + ip(p - l)wP-^{eZ)^. 

We denote the odd part and even part respectively by £^.42 and £^^2 with 

A2 = -A^foWx - {kl + kl)fQWx, 

A2 = {V'hfw.,^ - {kl + kl)xw^ + ^hWlhw^., + ]^p{p - \)wP-^{eZf. 

Then, we consider the problem 

-V2,xx + </52 - pw^" V2 = ^2 + -42 in R, 
1^2 (±00) =0, j tp2Wx dx = 0, 

as it is well known, which is uniquely solvable provided that 

(^^2 + A)w:r dx = 0. (3.57) 



(3.56) 



21 



(3.58) 



In fact, using the fact that w is an even function in the variable x, we have 

J •A.2'Wx dx = 0. 

On the other hand, we get 

/ A2WX d.T = -A^/o / wlAx- {kl + fc^)/o / wl Ax 
Jr Jr Jr 

= ^[A^f, + {kl + kl)fo] [ wldx. 

Jr. 

While ([338)) implies that 

J A2WX dx = 0, 

is equivalent to the following differential equation 

A^/o + (fc? + fc2)/o = onr. (3.59) 

In fact, we find that /o is identically zero by combining the equation (I3.59|) with the boundary 
condition ()4.1ip with the help of the non-degeneracy of F in ()1.3p . Whence, the solution to (|3.53p 
can be expressed as 

ip2ix,z) = ij2i{x,ez) +'4>22ix,ez), (3.60) 

where ip2i{x, sz) is in fact identically zero and -022 (a^, £z) is an even function in the variable x. The 
components in '>p2i{x, sz) and ^22{x, sz) are independent of the parameters /i and/2. 

In the same way, in order to cancel the error terms of order 0{s^) and improve the approxi- 
mation by solving problem p.54p . we collect all terms of order 0{s^) in 5'(V), which has the form 
£^^3 with 

A3 ^§3 + 83 + ^p{p - l)u.P-2(eZ)2 + p{p - l)wP-^{eZ + ^1)^2 + C3. 

We denote the odd part and even part respectively by and s^A3. As the arguments in solving 
p.56p . we need an orthogonality condition like p.57p . Hence, we compute the projection of ^3 
and A3 onto the kernel of the operator d^/dx^ — 1 + pwP^'^, which is spanned by w^- In fact, we 
obtain 

A3 = -A^hwx - (kf + kDhwx + C3 + 5(/o). 

Here g{fo) is a smooth bounded function independent of the unknown parameters /i and /2. 
Since the term A3 is even in the variable x and Wx is odd in the variable x, there holds 



A3WX dx = 0. 



R 
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On the other hand, using the same arguments as in p.58p . we get 

/ AsWxdx^—A^fi / dx — (fc^ + fcj)/! / Wxdx+ / C3Wxdx+ / gw^dx 

JR Jr Jr Jr Jr /oci\ 

(3.61) 

^-[A^h + {kl + kl)h] / wldx + diifo,e,ipi,ip2). 

Jr 

Here di is a smooth bounded function independent of the unknown parameters /i , /2 and the 
correction terms 1^93 and is Lipschitz continuous with respect to its parameters. By setting 

bi= I wl dx (3.62) 



we derive from (|3.6ip that 

/ AsWx dx = 0, 
Jr 

is equivalent to the following differential equation 



A^h + {kl + kl),h^^ onr. (3.63) 

Ol 

Now, we also choose /i by combining the equation p.63p with the boundary condition ()4.17p with 
the help of the non-degeneracy condition (|1.3p . The solution to p.54p can be expressed as 

ip3ix,z) = ip3i{x,ez) +ip32ix,ez), (3.64) 

where '4>3i{x,ez) is an odd function in the variable x and ■032(x,£z) is an even function in the 
variable x. The component in tp^i^x^ez) and '032(x,ez) are independent of the parameter /2. 



Remark 3.2. The recurrence procedure we described above is the same as the arguments in \35f . 
Whence, we can improve the local approximate solution to solve the first equation in \2.21\] up to 
0(e™) for any positive integer m. 



4 Boundary Correction Layers and Further Improvement of 
Approximations 

For any given (/2, e) G F, we have the inner expansion V in p.36|) . This approximation in general 
does not satisfy the boundary condition in (|2.21|) . In order to improve the approximation, we 
need to write problem (|2.2ip in local coordinates in the neighborhood of 9r e R"^, and then add 
boundary correction layers to the inner expansion p.36p . In fact, by recalling the scaled modified 
Fermi coordinates with variables 0, 77, s in Lemma 12.21 and also the translated variable x in p.33p . 
we will extend the local form of problem ()2.2ip to the infinite strip A in with notations 

A = { {x,9,ri) : X £R, (9, rj) e dB x R+}, 

(4.1) 

OA ^ { {x,0,r)) : X eR, e dB, r] ^0 }. 
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4.1 Local formulation of the scaled problem near boundary 

Here is the local form of problem (|2.21l) in the neighborhood of dT^ : 

Lemma 4.1. In terms of the coordinate system {s,6,r]), the equation in \2.21\) is expressed 



1 



1 



where B(u) = Bq{u) + Bi{u) + B2{u) with components 

C{e9) A{e6 



u + uP + B{u) = 0, 



Bo{u) 



. Mee) 



ll{ee)ll{ee) ll{e0)J " 



£^s(fc^ + k2)Us 



E{ee) R{e9)\ ^ . ^ 



(4.2) 



Biiu) 



I{e9) F{e6 



'h{ee)l2{ee) " Ii{e9)l2(ei 



-ug — 2es Usn 2es 



h{e9)l2{e0) 



h{e9)l2{e9) 



we have defined for e9 — i) 



570 (£0,0) 



> 0, 



AW 

cm 



92^0 (£0,0) 970 (£0,0) 



9p2 ' 


dp 


9^70 (£0,0) 


970 (£0,0) 


d-ddp ' 




9^70 (£0,0) 


570 (£0,0) 


dddp ' 


dp 



dp 

= ( 91 (£0,0), 



>o, 

970 (£0,0) 
dp 



m = 



i^(z9) = (gi(e0,O), 



9^70 (£0,0) 970 (£0,0) 
9t92 ' dd 

970 (£0,0) 



(4.3) 



dd 



and the differential operator B^{s,e9,er]) is of size 0{s*). Here, 01,02,61,62 are the coefficients of 
the operator B2 in h2. 28\) . 

The boundary condition in h2.21\) is recast as 



lesidv 97o(£0,O) \ ^ M(^e9) 



Z2(£0)2 \9p' 9p ' ^ /^(e0) 

_2 /(£0) /df 97o(£0,O)> 



q{s9)\dp' dp 



-^s'^Us + e'^bi{e9)s'^Urj 



(4.4) 



where we have denoted 



M{e9) 



+ e^65(£0)s^Ms + £266(£S,£0)Me + D{u) = 0, 

970 (£0,0) dv\ /d-fo{e9,0) dv 



dd 'dpi \ dp ' dd I 
The differential operator D{es, e9) is of size 0{e*) and functions bii{e9), b^{e9), bQ{es,e9) satisfy: 

|64(£0), 65(£0), 66(£s,£0)| < C{l + \es\''). 
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Proof. The reader can refer to Section 4 in [34 and the references therein for the details of the 
proof. □ 



As we have done in p.33p . we introduce a further changing of variables and define a new 
function v{x,9,ri) as follows 

2 



i{s,9,r]) = v{x,9,r]) with x — s — ''^^ fi{ed , erj) , 9 = 9, ?? = Jy- 



(4.5) 



We now want to express the problem in the new coordinates. Whence we need the following 
formulas 



2 

Use = [ - 



I j ^xx ~r Vxrj , 



uee 



2 2 

2 = 1 = 

2 2 ^ 2 

2=0 1=0 'i=0 



i=0 j=0 j=0 

Using Lemma 14.11 we get the local form of the problem (|2.2ip in the new coordinates and then 
extend it to the infinite strip A 



S{v) = + ,2} n\ '"ee + wr^^vv + + B{v) =0 in A, 



where we have denoted 



q{e9y'"' qis6 



B{v) =Bo{v)+Bi{v)+B2{v), 



(4.6) 



(4.7) 



liis9) 



j=0 



j=0 



i=0 



A(e6i) 



Z^(ei 



2 

E^'^'^p/*) + ( - E^'^'^p/^)^- + ( - E^'^'^p/' 



C{e9) A{e6 



Il{e9)q{e0) ll{e9) 

2^2 2 



1=0 



4=0 



i=0 



£;(e6') i?(£6 



^?(£0)^i(£^^) il{s9) 



z 



i=0 



e^b{e9,er])ix + '^e'fi U^, 



i=0 
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'hieOMee) 



i=0 



F{e6 



z 



■ Vg 



4=0 



4=0 

2 



i=0 ^ 



i=0 
2 



(e0)/2(£^) 



j=0 













1- We 




[( 



j=0 



4=0 



4=0 



2 



The boundary condition is 
2 



4 = 



2£ /diy 970 (£61,0) 



/(£0) /diy d-fQie0,O) 



4=0 



4=0 



ll{ee)\dp' dp 



[x + J2 £74) + e%{ee) [x + J2 £^4) 



4=0 



4=0 



4=0 4=0 



(4.8) 



2 

+ £2&6(X, £0) [( - E £'^'V'5/«)«- 

4 = 

4.2 Boundary correction layers 



ve 



+ D{v) = on OA. 



In this subsection, we will find some boundary correction layer terms, say i^i's, by also the recur- 
rence method. The method is basically the same as that in [5S] and [TO] . 
By recalling p.36p . we take 

3 

ui =V = w{x) + seZ + Ye'-Lpi, (4.9) 
4=2 



as the first approximate solution of the problem (|4.6p and (|4.8p on A. Then we compute 

S{ui) ^ e^A^eZ + sXaeZ + B3{ui) in A. (4.10) 
On the boundary, the errors become 

$ = £$1 + £^$2 + £^$3 + D{ui) =0 on OA, 
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where we have denoted 

$1 = (^Iie9)x + I{e0)fo + Vpfo)w,, 



I{ee) /du djo{e9,0) \f 



2 M{e9) 







M(£6') 



(a^ + /o) (V 



To improve the approximation, we need to cancel the terms $i's by adding boundary correction 
terms to the inner approximate expansion in (I4.9p . This can be done step by step. On the boundary 
it is natural to take 



Vp/o + I{ee)fo = 0, 



(4.11) 



which will lead to the cancelation of {I{e9)fo + pfo)wx in the components of the error of first 
order of e. Recall that by combining of p.59p and (14. lip we have chosen /q = in Section [3] On 
the other hand, to cancel the first order term like sl{s9)xwx on the boundary, we shall introduce 
a boundary layer term, say (jji. Indeed, we first introduce the term 



<j)ix,9,7j) = eAoie9)b{VXol2{s9)7j)Zix), 

with (j){x, 9, rj) satisfying 

ixx + ,2} n\ ^se + ,2} a\ ^m -4> + pwP^^4> = 0{e-^) in A, 



drj 



-Co{s9)Z, on dA, 



where co{e9) is a function in the parameter 9 of the form 

CQ{e9) — I{e9) / xWxZ dx. 
Jr 



(4.12) 
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In fact, we can choose 

4>{x,d,f^) - e sm{^ol2{e9)v)Z{x) = e^ii{x,d,,j). (4.13) 

Then, by CoroUarv lA.l.ll there exists a unique solution (denoted by 0i2) of the following problem 

1 I 

+ ,2/ Q-. 4>12,'q'q " 012 + PW^" Vl2 =0 in A, 



= I{e9)xwx — co{e9)Z on dA. 



dr] 

Moreover, (j)i2 is even in the variable x. By choosing a smooth cut-off function xo in the form 

Xo(i) = l if |i|<l, Xo(t)=0 if \t\>2, (4.14) 
we can set the first boundary layer term by 

(1)1 ^ £Xa{£v)'^i{x,0,r]) with ^i{x,e,f]) ^ (f>ii{x,9,r]) + (j)i2{x,e,'i]). (4.15) 

Thus we finish the first step. 

To proceed the second step, let U2 = ui + 0i be the second approximate solution. We again 
compute the new error 

Siu2) - Sim) +2e' *M + ^'ifP) ^0 *i + ^^^j^w^-'ieeZ + 



C{ee) A{ee) I{ee) 



ll{ee)ll{e6) l\{ee) h{ee)h{ee) ) ^^^^q^ 



XO*l,r, 



2/1 ^ ^M) \ X 

+ Co(0ii) + Lo(0i) + A^o(0i) in A, 
where Co{(j)ii), io(0i) and No{(j)i) are of size 0{e^) and S'(ui) is defined in (|4.10|) . 

Co(0n)-5(0n)+e (^^^^-^^ - ^ - ^-^-^^ j 



+ ^\lIP)+'"^l(i2(s^)j^"*^-'" 



Cjed) 

1 ^ ^(e^*) 



Lo(0i) = S(0i) - £ ( ^ipyiipy - - /,(e0)^,(e^),)^°*i- 
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Note that /o = 0. On the boundary, the new error becomes 

(^pfi + I{ee)h + I{ee)e xZ^ - Vpe z \ + e'l{eB)x </>i2,: 



q{e9)\dp' dp 



iliee) 



'12,6 



Vp/2Wx + I{£0)f2Wx + Vp/i(eZj; + (t>i2,x) + I{e9)fi{eZx + cj)i2,x) ~ V2,-n 



qie0)\dp' dp 



x'^{eZx + (f'i2.x) + 2fixwx 



+ e^b4{ed)x^'i'i^,^ + e^b^{£e)x^Wx - £^h&{x, ed)Wi)foWx 
-£^^^^x(V^/iu;x- V^eZ) +5(u2) = on 9A. 

In order to cancel the error terms of order 0(£^) on the boundary, we do the same thing as that 
in the first step. We first set 



Vp/i + I{e9)h = 0, 



(4.17) 



and 



Then we choose 



JR L 



xZx — V„e Z 



Zdx ^ Wpe+-I{ee)e = 0. 



(4.18) 



'Ci{e9) ~ C2{e6 



where ci (eO) and C2 {s9) are functions in the parameter 9 of the form 

M{e9) 



m{^\nl2{e9)rj)Z{x) = e^cj)2iix,9,r]), (4.19) 



ci{e9) = I{e9) / x4>i2.xZdx, C2(e6 
Jm 



■ gZdx. 



(4.20) 



Now, by Lemma Fa. 1.31 there exists a unique solution (denoted by (/)22) of the following problem 



ll{e9) 



022,6 



q{s9) 



i>22, ip-, - (t>22 + pw^ ^(f>22 = hi in A, 



90; 



'22 



dri 



= I{e9)e xZx -WpeZ + I{e9)x(f>i2,x - ci {e9)Z 



where hi is a function in the form 
1 



q{e9)' 



hi =-2 ^^x'o ^i,v ~ ^^^j^wP-^{eeZ + cj^^f 



C{e9) A{e9) 



1(69) 



Il{e9)ll{e9) q{e9) h{e9)l2{e9) 



M{e9) 



'12,9 - C2{e9)Z on 9A, 



2x 



I{e9) 



q{e9) h{e9)l2{e9) 
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Moreover 022 is even in the variable x. By recalling the cut-off function in (|4.14p . we set the 
second boundary layer term by 

02 =e\oN)*2(2:,f,?7) with '^■2ix,9,r]) ^(t>2i{x,e,r])+(l>22ix,9,r]). (4.21) 

For the completeness, we here give the third step although it is the same as we have done in 
the above. Let U3 = U2 + 4'2 be the second approximate solution. We compute the new error 

^(.3) =^(.2) + ^-^^X^'*! + 2 j±^^x'o *2„ + e^^x^' ^2 

/ 1^ P-2 sf C{ee) A{ee) I{ee) \ ^ 



- e3x(fc? + A:|)xo*i,^ +Co + Ci + Li((/)2) + A^i(02) in A, 
where Ci((/)2i), io(02) and No{(j)2) are of size 0(£^) and S'(u2) is defined in (|4.16p . 



r. 3 Fi^S) . 



h{ee)l2{ee) 



C{ed) A{ee) I{e9) 



iVi((/)2) =(U2 + (p2y -ul- pwf- V2 - V{V - \)w^-^eeZ^2- 
On the boundary, the new error becomes 

+ e^/(e6')a; 022, + e'h^{eQ)x^'^x,r) + e^65(e6')a;^u;a; 

x(Vi;/iw^ - VijeZ) + 5(7/3) =0 on dK. 



q{e9) 

In order to cancel the error terms of order O(e^) on the boundary, we first set 

Vp/2 + ne9)f2 + 9\{9) = 0, (4.22) 

by introducing the term 

9l{9) = ^T^V^e [ xZw^dx / [ wldx. (4.23) 
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H^,e,rj) ^ e^^^f-^8m(^h{e0)T])z{x)=ec^3i{x,e,T]), (4.24) 
\/Xnh eO) V / 



Then we also introduce the term of the form 

where c^^eO) is a function of the form 

csiee)^ I C4ix,ee)Zdx. (4.25) 

In the above, we have denoted 

C4(x, eO) ^I{e9)x ip2,x + I{e9)x 022,1 - <(52,r, + hi{e9)x^ 
+ h{ee)x^ - ) ' Vahxw^, 



(4.26) 



Again, by Lemma [A. 1 .31 there exists a unique solution (denoted by ^32) of the following problem 



C4{x,e9) — C3{s9)Z on DA, 



drj 

where /i2 is a function of the form 



, C{e9) A{e9) I{e9, . ^ , p_2 



+ Co/e^ + x{kl + kl)xQ'ifi,x- 
Moreover (/)22 is even in the variable x. We can set the third boundary layer term by 

'?!>3 =£^Xo(e?7)^'3(a;,6',77) with *3(x, 6*, 77) = <?!)3i(a;, 6*, r?) + ^32(2;, 6*, ??)• (4-27) 

Remark 4.2. In fact, we can proceed the above arguments step by step to find boundary correction 
terms and get rid of the error terms up to order 0(6™) for any positive integer m. 

4.3 Summary 

We conclude that for any given parameter pair (/2,e) S F, our final approximate solution to 
the problem (I2.2ip near the surface is expressed in the local form by 

3 

U4 = w{x) + eeZ + ^ + 01 + 02 + 03- (4.28) 

For a suitable perturbation term 0, if we locally set U4 + as the solution to problem (I2.2ip . 
the problem can be recast as follows 

S{U4 + 0) = S{u4) + L2(0) + B(0) + iV2(0) = 0, (4.29) 
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where the linear and nonlinear operators are in the form 

N2{(i)) = + <i>r -u\- V, 

with boundary condition 

5o(0) - 0^ + 5(u4 + 0) =g. 

The error of the approximation is 

Ey =5(U4) 

= i'i^eZ + eMeZ + ( B4 + B4 + ^4 + 64 + ©4 + ©4 ) 

6 6 10 

i=4 s=4 i=5 

In the above, we also have denoted 



DQ{(t)) =el{ee)x(j)^ + g ,2, ' x(t)0 + e'^h^{ee)x-^ (t)^ 



eO) + V pfi ) fi+2,z 

2 2 ^ 

i=l i=2 

2 2 2 4 



j=l i=l i=2 

2 2 4 3 



i=l i=l 1=2 j=2 

1^ i=l i=2 i=2 



4 3 



+ e%{x, sO) (sWi^eZ + E ^'Vi,e + E ^Vi2,0) 

i=2 i=l 
2 4 3 



i=l i=2 i=l 

For later use, we decompose the error as two components 

El = Ell + E12, 

where we have denoted 

£11 = e^A^eZ + eXoeZ and £^12 = Ei - En. 
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4.4 Size of the errors in weighted Sobolev norms 

To estimate the size of error, we have to introduce some suitable weighted Sobolev norms. Here 
we use the same norms as those in 12,: for a function h{x,z) defined on a set E G R"^, and for 
< Q < jig and 4 < q < +oo, we set 

||^||<;,e;E= sup e^l^'l||/l||i5(B((a;^2)^i)), 

(x.z)eE 

2 (4.36) 

Here B{{x, z), 1) denotes the ball of radius 1 centered at (x, z). 

For the application of contraction mapping theorem in the procedure of finding the perturbation 
term 0, we need analyze the properties of En and E12. The reader can refer to Section[71 From 
the uniform bound of e in p.37p , it is easy to see that 

\\En\\g,,.,e<Cei+'-l. (4.37) 

where & is defined in p.29p . 

All terms in £'12 carry in front, we then claim that 

Il^^l2||g,e;e <C£5+4-f. (4.38) 

A rather delicate term in E12 is the one carrying A'"/2 since we only assume a uniform bound on 
II /2||L9(r)- For example, we have a term Ki = f2Wx in S{w) which has bound like 

||ifi||g,,;e<Ce^+4-f. 

Similarly, we have the following estimates 

||5ll,,,;a6 <Ce^+'"'. (4.39) 

Other terms can be estimated in the similar way. Moreover, for the Lipschitz dependence of the 
term of error £^12 on the parameters /2 and e for the norms defined in p.32p and p.37p , we have 
the validity of the estimate 

||i^i2(/2,e)-£;i2(/2,g)||,,e;e < Ce^+^-f ( II/2 - + ||e - g||fc )■ (4-40) 
Similarly we obtain 

Il5(/2,e)-g(/2,e~)||,,,;ae < Ce^+^-f ( II/2 - /zlU + ||e - gHb)- (4.41) 
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5 The Gluing Procedure 

Recall that, in Sections [3] and HI we consider problem (j2.2ip in a small neighborhood of and find 
a local approximate solution. In this section, to get a real solution to (|2.2II) by the perturbation 
method, we use a gluing technique (as in [10^ ) to reduce the problem in D,^ to a projected problem 
on the infinite strip 6(cf. ((3?29| ) in R^. 

Let a < ro/100 be a fixed number, where tq is a constant defined in (|2.7|) . We consider a 
smooth cut-off function ricr{t) where t e IR+ such that ricr{t) ^ 1 ioi Q < t < a and ri{t) ~ for 
t > 2(7. Set Ti^{s) = ?yCT(£|s|)i where s is the normal coordinate to F^. Let 1*4(5,2) denote the 
approximate solution constructed near the surface F^ in the coordinates (s, z). We define our first 
global approximation to be simply 

W = vlAs)u^- (5.1) 

Obviously, is a function defined on Q^, which is extended globally as outside the 6a /e- 
neighborhood of F^. 

For u = W + (j) where (j> globally defined in fig , denote 

S{u) = AyU — u + in rjg. 

Then u satisfies (|I.ip if and only if 

£(0) = ~E~ N{4)) in fie, (5.2) 



with boundary condition 



where 



dW 

on dn„ (5.3) 



N{4>) = {w + 4>T -wp- pwp-^cj>. 

We will look for (f> in the following form 



where, in the coordinates (x, z) of the form p.33p . we assume that (p is defined in the whole strip 
&. Obviously, (|5.2p - (|5.3p is equivalent to the following problem 

r,ls [Ay^ - + pWP''^) = r^l [-Niv'^g^ + ^j) - E - pM^f" V] , (5.4) 



34 



(5.5) 



On the boundary, we get 



|^ + (l-.l)|^+^^^-0. (5.7) 

The key observation is that, after solving (|5.5p and (15. 7p . the problem can be transformed to the 
following nonlinear problem involving the parameter tp 

C{(l>)=r,l\ -N{r,ls(l) + i;)-E~pWP-'ij] in e, (5.8) 



dd) dW 

7r^+vljr^=0 on 96. (5.9) 
arie oris 

Notice that the operator C in fi^ may be taken as any compatible extension outside the 6a /e- 
neighborhood of T/e in the strip 6 and the operator d/dn^ may be taken as any compatible 
extension outside the 6cr/e- neighborhood of T/e on the boundary d&. 

First, we solve, given a small <j), problem (|5.5p and (|5.7I) for tp. Assume now that (j> satisfies the 
following decay property 

|V(/.(y)| + |0(y)| <e-^/^ if \s\>a/e, (5.10) 

for certain constant 7 > 0. The solvability can be done in the following way: let us observe that 
W is exponentially small for \s\ > cr/s, where s is the normal coordinate to T/s. Then the problem 

AV' - [1 - (1 - riDpWP-^] ^ = h in 
dip dW drils 

7i— = -i-%)t^ ^ir^<P onar^e, 

has a unique bounded solution Tp whenever ||/i||oo ^ +00. Moreover, 

<C\\h\\^. 



Since N is power-like with power greater than one, a direct application of contraction mapping 
principle yields that (|5.5I) and (|5.7I) has a unique (small) solution ip = "ipi^i) with 

m4>)\\L^ < Ce[ ||0|U~(|«|>./e) + l|V0|U»(|,i>^/,) +e--/^ ], (5.11) 
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where \s\ > a/e denotes the complement in fig of cr/e-neighborhood of F/e. Moreover, the nonhnear 
operator tp satisfies a Lipschitz condition of the form 

UiM - iji(j}2)\\L^ < Ce[ \\(bl-h\\L^i\s\>a/e} + 1 1 V(/.i - V02 1 1 L~ (|s|>a/e) ]• (5.12) 

Therefore, from the above discussion, the fuh problem has been reduced to solving the following 
(nonlocal) problem in the infinite strip & 

C2{<P)^vl[-N{7^ls<t> + ^m-E-pWP-'i^{(t>)\ in©, (5.13) 
dW 

'^{<t>)+Vs^ = onde, (5.14) 
one 

for (p e W'^''^{&) satisfying condition (jS.lOp . Here £2 denotes a linear operator that coincides with 
£ on the region \s\ < Scr/s, *B denotes the inward normal derivatives of & that coincides with 
inward normal d/dn^ of fl^ on the region \s\ < Sa/e. 

The definitions of these operators can be shown as follows. The operator £ for \s\ < Safe is 
given in coordinates {x,z) by formula (|3.33p . We extend it for functions defined in the strip & 
in terms of (a;, z) as the following 

£2(0) =i2(0)+x(e|a;|)B(0) in6 (5.15) 

where x(^) is a smooth cut-off function which equals 1 for < r < 10a and vanishes identically 
for r > 20(7, L2 and B are the operators defined in (|4.30p and p.34p . Similarly, the boundary 
conditions can be written as 

xis\x\)Do{<P)-^+x{£\x\)DiW + (t>)=x{£\x\) 9 on 96 (5.16) 

where the operators Dq and D are defined in (|4.3ip and (|4.8p , is the unit inward normal of d& . 

Rather than solving problem ()5.13|) - (|5.14|) satisfying the boundary condition, we deal with 
the following projected problem: for each pair of parameters /2 and e in F, finding functions 
(j) e W^'i{e), c,de i«(r) and Ai, A2 such that 

£2(0) = -xEi - X^^3('^) + c{ez) xwx + d{ez) xZ in 6, (5.17) 



xDoi^) - JT^ + XD{W + cP) ^ X9 on 96, (5.18) 

•{x, z)wxdx — Ai, / (f> {x, z)Z{x) dx — A2 in Fe, (5.19) 
Jr 

where N3{(j)) = N{ri^g<j) + tp{(j)))+pWP~^ip{(l)). In Sections [6] and [3 we will prove that this problem 
has a unique solution (f> whose norm is controlled by the |1 • ||,j.g-norm, not of the error component 
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Ell, but rather of the components £'12 and g. Moreover, (j> will satisfies (|5.f Op . The reader can 
refer to the conclusion in Proposition 17. II 

After this has been done, our task is to adjust the parameters /2 and e such that the functions 
c and d are identically zero. It is equivalent to solving a nonlocal, nonlinear coupled second order 
system of differential equations for the pair (/2,e) with suitable boundary conditions. Indeed, we 
will derive the system of differential equations for the unknown functions /2, e in Section [SJ and 
then show the solvability of this system on the infinite dimensional space F(cf. p.38[) ) in Section 

m 



6 The Invertibility of £2 

Let £2 be the operator defined by (|5.I5p and g be the functions in (|4.34p . Note that the function 
is even in the definition of £2- In this section, We study the following linear problem: for 
given h € L9(6) and g e Li{d&), finding functions <j> e W'^-'^{&), c,de Li{dG) and Ai,A2 such 
that 

^2i(t>) = xh + c{ez)xwx + d{ez)xZ in 6, 
xDoi^) ~ ^ + xDiW + qj) ^ X9 ond6, (6.1) 
/ (f>Wx{x)dx — Ai, / (f>Z{x)dx — A2 in T^. 

Proposition 6.1. // a in the definition of C2 is chosen small enough and h G L'^{&), then there 
exists a constant C > 0, independent of e, such that for all small e, the problem id.l]) has a unique 
solution (j) = Ti(h,g) with suitable Ai and A2 which satisfy 



(6.2) 



< C(||/i||l,(6) + ||g||L5(ae)), 
||Ai||L<,(r,) < C'(||/i||i,(e) + ||g||L9(ae)), V z = 1,2. 
Moreover, ifh,g have compact supports contained in \x\ < 20a /e, then 

\(l>{x,z)\ + |V^(a;,z)| < ||(/)||loo e-^"^" for \x\ > AOa/e. (6.3) 

Proof. Note that the problem can be written as 

<Px. + A^^c/) - + V - - p{W^'-^ ^wP-^)cj)- xB(cj>) 

+ xh + c{ez) xwx + d{ez) xZ in 6, 

d<l) 



. -X9 + xDo{<l^) + XD{W + (j)) ond& 

I (pWxdx = Ai, / (f>Z{x)dx — A2 in T^. 
Jr Jr 
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Let 

^^T,(^xh-p{WP-'-wP-'^~xBW, g), 

where 

G(0) = -X9 + x5o(0) + XD{W + 0), 

and Ti is the bounded operator defined by Lemma [6.2.41 
The key point is that the operator 

^4(0) = -xB{cb) - p{WP-^ - wP-^)cj,, 

is smah in the sense that 

||B4(<^)||,,,< Gallon,,,. 

Similar results hold for G(</>). Hence, the results can be derived by the invertibility conclusion of 
Lemma [6.2.41 if we choose a sufficiently small. 

Since x is supported on \x\ < 20a/e, then (j) satisfies for > iOa/e a problem of the form 

+ Ar'0-(l + o(l))0 = for \x\>40a/e, ^ ^ ^■ 

Hence, the validity of formula (|6.3p can be showed easily. □ 



7 Solving the Nonlinear Projection Problem 

In this section, we will solve (|5.17p - (|5.19p in (3. A first elementary, but crucial observation is the 
following: The term 

Ell e^A^eZ + eXoeZ, 

in the decomposition of Ei, has precisely the form d{ez)Z and can be absorbed in that term 
xd{ez)Z . Then, the equivalent equation of (|5.17p is 

^^2(0) = xEi2 + xNz{4') + c(ez) x^x + d{ez) xZ. 

Let Ti be the bounded operator defined by Proposition 16. II Then the problem (j5.17p - (j5.19p is 
equivalent to the following fixed point problem 

= Ti{xEi2+xNm, Xg) = -4(0). (7.4) 

We collect some useful facts to find the domain of the operator A such that A becomes a contraction 
mapping. 
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The big difference between En and E12 is their sizes. From (|4.37p and (|4.38l) 



\E- 



while £'11 is only of size 0(e^+^/^ •^/*). Similarly, we have 



(7.5) 



(7.6) 



The operator Ti has a useful property: assume h has a support contained in < 20a/e, then 
= Ti{h) satisfies the estimate 



\(f>\ + |V0| < M\oo,s e"'"/" for > 40a/e. 
Recall that the operator tp{(l)) satisfies, as seen directly from its definition 



(7.7) 



< Ce 



I'/'I + IV./.I 



lL°°(|2;|>20cr/£) 



and a Lipschitz condition of the form 



\L°°{\x\>20<T/e) 



Now, the facts above will allow us to construct a region where contraction mapping principle 
applies and then solve the problem (|5.17l) - (|5.19l) . Consider the following closed, bounded subset 



e i?2(6) 



101 + |V0| 



i+4- 



< 



(7.8) 



L°°{\x\>iaa-/e) 

We claim that if the constant c, is sufficiently large, then the map A defined in (|7.4p is a 
contraction form 2D into itself. Let us analyze the Lipschitz character of the nonlinear operator 
involved in A for functions in 2) 

xiV3(0) = xNii^ + V'(0)) + xpW^p- V(0) 
Note that Ni{(p) = p[{W + t(p)P~^ - WP-^]ip for some t G (0, 1). From here it follows that 

iiVi(^)i<c(i(^r + i^n. 

Denoting S'o- = © n {|a;| < lOcr/e}, we have that for ^ e D 

imm,,, < c[ \m^,,^, + \m%^, + \\m\\i,,-^s. 

Using Sobolev's embedding, we derive 

<C{ 



2 



P 

qp,Q 



2 

2q,0 



P 



2 

2,q,0 
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Using estimates (|5.1ip . the facts that e £>, (I7.7p . that of the area of Sa is of order 0{cr/e) and 
Sobolev's embedding, we get 



p 



2 < r'p-'^/4e 



2 

2,9, e 



Hence, from the properties of W and ^{(j)) we obtain 

As for Lipschitz condition, after a direct calculation we find 



(7.9) 



'2j|k,e 



< C 



WlWqp^e + ll'PlIhg.e + ||V'2||^pj + ||¥'2||29,e 



X (||<Pl - ^2\\qp,e + IIV'l - ¥'2||29,e)- 



Hence, 



\\N^{4>) - N^m\q,9 < 11^1 (0 + ^(0)) - + ^{mWe:.s. 

< - 4>\Wg;S. + II'/' - 0ll<;p,0;S„) 

+ w(||V(0) - V'('^)I|29,,;S. + ||V'('^) - i'mUp,,;S.), 

where v = vi + V2 with 

= \\Mfq;'s;S. + IIV'C'^OII^p^U + ll'^'l|29,,;S„ + ||V'(^OI|29,e;S„. 

Arguing as above and using the Lipschitz dependence of '0 on (j), it can be derived 
||xiV3(0) - xN^mq., < C(£(H4-f + „ 



(7.10) 



Now, we can find the solution of (|7.4p in the sequel. Let G 2) and = .4(0), then from 
(1731)- (d]) and dZSD 



rll2,9,e 



< llTi 



(|+4-f)P I ^^^2^9- 



Choosing any number <^ > C*||Ti||, we get that for small e 



J^\\2,q.Q < 



+4- 



From (O 



kl + |V^| 



< lli^lloo < \W\ 



2,9,e e 



L°°{\x\>40cr/e) 

Therefore, v £ Tl. A is clearly a contraction thanks to (j7.10p and we can conclude that (j7.4p has 
a unique solution in J). 
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The error £'12 and the operator Ti itself carry the functions /2 and e as parameters. For future 
reference, we should consider their Lipschitz dependence on these parameters. (|4.40[) is just the 
formula about the Lipschitz dependence of error E12 on these two parameters. The other task can 
be realized by careful and direct computations of all terms involved in the differential operator 
which will show this dependence is indeed Lipschitz with respect to the M^^'''-norm (for all e). 

Within the operator, consider for instance the following term involving f2 

Then we have 



IIQ/.(0)IIL(B)<^''"' / |A^/2r(sup / \Mx,zWdx 
Let ^{z) = z)\'^dx. Then there holds 

supAi(z)<£ / |0,|«dx + i / |0,|«-i|V^0,|dx 

< isupAi(z) + -^ / iV^^.I'dx, 



(7.11) 



and we can obtain 
Therefore, 

||Q/.(0)||,,,<C£||/2|U. 

Similar estimates can be applied to other terms in the operator involving A'"/2. 
For the linear operator Ti , we have the following Lipschitz dependence 

\\Tlif2)-T,ih)\\2^,^S<Ce\\f2~f2\\a- 

Moreover, the operator A^3 also has Lipschitz dependence on (/2,e). It is easily checked that for 
G 2? we have, with obvious notation 

||xA^3,(/„e)(0) - X^3,(/2,e)Wll9,. < CeH^-f [ \\f^ _ f^\\^ + ||e _ 

Hence, from the fixed point characterization we get that 

Mf2,e) - Hf2,e)\\2,g,e < Ce^+^-f [ II/2 - hWa + \\e - e\\, ]. (7.12) 

As a conclusion, we give the proposition. 

Proposition 7.1. There is a number <7 > such that for all e small enough and all parameters 
(/2,e) in F , problem \5.n\j - l5.19\) has a unique solution (p ~ 0(/2,e) which satisfies 
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101 + |V0| 



< 



2,9, 



-(J je 



L°=(|a;|>40cr/e) 

Moreover, (f)(f2,e) depends Lip schitz- continuously on the parameters md e in the sense of the 
estimate (TJ^. 



□ 



8 Estimates of the Projections Against Wx and Z 



As we mentioned in Section [51 we will set up the system of differential equations for the unknown 
parameters /2 and e defined on T which arc equivalent to making c, d zero in the system (jS.f 7p - 
(|5.19[) . On the boundary, we have imposed the conditions (j4.22p and (I4.18P for the parameters /2 
and e, which are restated respectively in the following 



Vp/2 + I{ee)f2 + m{e) ^ O, V^e + -Iie9)e = 0. 



(8.1) 



On the interior of F, these equations are obtained by simply integrating the equations (j5.17p (only 
in x) against Wx and Z respectively. It is easy to derive the following equations 



+ xNsicj)) +A^'(j> + xB(0) + p{WP-^ 



Wxdx — 0, 



xEi + xN3{^) + A^^0 + Ao0 + xBW + p{WP-' - w^-'] 



Zdx = 0. 



(8.2) 



.3) 



where the error term Ei is defined in (|4.32l) . the operators and B are defined in (|5.17p and 
()3.34p . It is crucial to estimate the terms 



EiWxdx and / EiZdx 



The same arguments can be applied to other terms in (j8.2p and 
for the components in (|8.2p and (|8.3p into several parts. 



Now, we divide the estimates 



8.1 Estimates of the projection against 



First, multiplying (j4.32l) by Wx and integrating over the variable x, using the decomposition of Ei 
in (|4.35p and the fact that Wx is an odd function in x, we obtain 



EiWxdx 



Ei2Wxdx. 



More precisely, there holds 



/ Ei2Wxdx^s'^ / [B4 + G4 + S4 + pip - l)w''^^eZ(j)3]wxdx 

» 6 6 10 



Wxdx. 
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By direct calculation, we derive that 



10 (8.4) 

1=5 



where 71 , 72 are two constants and 

cTi = / wldx. 



Here and below we denote by = i,'2,i, generic, uniformly bounded continuous functions of 
the form 

he = bUf2,e,V^f2,V^e), 
where additionally bi^ is uniformly Lipschitz in its arguments. 

8.2 Projection of terms involving 

We will estimate other terms that involve (j) in (18. 2p . 

xN2{(p) + A^^<P + xB{(l)) + p{WP-^ ~ wP-^)(l)'\w^dx. (8.5) 

Using the condition in (15. 19^ . we first estimate J-g^A^'(j)'Wxdx, the estimate for the term can be 
done as follows 

Ti(z) = / A^'(j)w^dx ^ A^' / (j)w^dx = A^^Ai. 
Jr Jr 

With the help of the Proposition 17. 11 we have the following estimate 

l|Ti|U,(r,)<Ce^+^-t. 

The last two components in (|8.5p are 

T2 = / B{(l3)w^dx and T3 / p{W''-^ - wP-^)(l)w^dx. 
JR Jr 

Here we recalled the definitions of the operator B in (13.341) and the local approximation W in 

(|5.ip . We make the following observation: all terms in B{(l>) carry and involve powers of x times 

derivatives of 1, 2 or two orders of 0. The conclusion is that since has exponential decay then 



Hence there holds 



l^\T2ieWd9<Ce^m,,,,- 
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In B{(f>) we single out two less regular terms. The one whose coefficient depends on A'"/2 explicitly 
has the form 

•'^ [ 1=0 } X 

Since (j) has Lipschitz dependence on (/2, e) in the form (|7.12p . we see that 

||T2*(/2,e)-T2*(/2,g)||L.(r) < (II/2 - /2||a + ||e - elb). (8.6) 

The other arising from second derivative in y for is 



1- (l + £(x-^e'/i 



-2 

I 



1=0 



dx. 



T2** = / A (j)'Wx 
Jm 

We readily see that 

llT2,,(/2,e)-T2,,(/2,g)|lL,(r) < Ce^+^+l (y/^ - /2IU + ||e - glU). (8.7) 

The remainder T2 — T2* — T2** actually defines for fixed e a compact operator of the pair (/2, e) 
into L'^(r). This is a consequence of the fact that weak convergence in W'^''^{&) implies local strong 
convergence in W^''^(&). If f2,j and Cj are weakly convergent sequences in W'^''^{&) then clearly 
the functions (f'{f2,j,^j) constitute a bounded sequence in W^''^{&). In the above remainder one 
can integrate by parts if necessary once in x. Averaging against Wx which decays exponentially 
localizes the situation and the desired result follows. 
Let us consider now the term 

T3(z)= / p{WP''^ ^wP-^)(t)Wxdx. 
Jr 

3 3 

Since the term W ~ w{x) + eeZ + ^Vj + 't'i can be estimated as 

i—2 2—1 

3 3 

e\e{ez)Z{x)\ + kV.I + E I'^^l ^ ^^(^ + 1^1')^"'''. 

1=2 1=1 

we easily see that for some k > the uniform bound holds 

From here we readily find that 

||T3|U,(r) < Cet||0||2,,,e < Ce^+^+i 
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We observe also that the term in (jS.Sp such as 

can be estimated similarly. In fact, using the definition of N^ic/)) and the exponential decay of Wx 
we obtain 

\\r4LHr)<CMi,^,<Ce'-t 
These terms define compact operators similarly as before. 

8.3 Estimates of the projection against Z 

We observe that exactly the same estimates can be carried out in the terms obtained from inte- 
gration against Z . So the remaining thing is to compute the term EiZdx. 

Multiplying ()4.32p by Z and integrating over the variable x and using the decomposition of Ei 
in (|4.35p . we get 

J EiZAx^ j EiiZdx + J E12ZAX, 

where 

/ EiiZAx = e(e2A^e + Aoe) / Z^dx = e^A^e + eAoe. 



On the other hand, we have 



( Ei2Zdx = e'^ [ [Bi + Gi+t>i +p{p-l)wP'^eZ(f>3]Zdx 

» 6 6 10 

+ / [Y,£"S^ + Y,£"T, + Y,e'G^+/'B5\zdx. 

jR 



The components in B4,G4and2)4 are even functions in the variable x and independent of the 
parameters /2. Here, the function 03 has the form 

(/)3(a;, z) = ip3i{x,ez) + ^'32(2;, ez). 

The components in ip^i and '032 aie independent of the parameters /2. Moreover, -031 is an odd 
function in the variable x and ^^32 is an even function in the variable x. Therefore, adding up all 
terms together, we conclude that 

. 10 

/ E,2Zdx ^e^lA^e + e^blA^ f2 + y^e%,if2,e,\/^ f2,S/^e). 

Here we denote by b\^,j — 1,2 and bi^, generic, uniformly bounded continuous functions, moreover, 
hie is uniformly Lipschitz in its arguments. 
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9 Solving the System for (/2,e): Proof of Theorem HJJ 
9.1 Proof of Theorem [TT] 

Using the estimates in previous section, we find the following nonlinear, nonlocal system of differ- 
ential equations for the parameters (/2, e) in the variable y = ez G F 

Cl{f2) = A^f2 + {kl + kj)f2^jie + j2e^A^e + A'h, in T, (9.1) 
C2,{e) = -e^A^e-Xoe^M2e in T, (9.2) 



with the boundary conditions 



^ + //2 + $n = on ar, (9.3) 

OT 

1^ + i/e = on ar, (9.4) 
OT 2 



where 71 and 72 are two constants defined in (|8.4p . / and D\ are smooth functions defined in (|4.3p 
and (|4.23p and t denotes the inward normal of dT. The operators Mie and M2e can be decomposed 
in the following form 

Mu{f2,e) = AUf2,e) + Ki,{f2,e), 1 = 1,2 

where Ki^ is uniformly bounded in L'^{T) for (/2, e) in F and is also compact. The operator is 
Lipschitz in this region, see (j8.6p - (l8.7p . 

P/e(/2,e)-Ai,(/2,g)|U,(r) <Ce5+i[ 11/2-/211, + lie- (9.5) 

Before solving (|9.ip - (|9.4p . some basic facts about the invertibility of corresponding linear oper- 
ators are in order. We first consider the following problem 

£I(/2) = A^/^ + {kl + kl)f2 = in r, ^ + 7/2 = on dV. (9.6) 

OT 

Lemma 9.1. Under the non- degeneracy condition 0/ T in if h Cz L'^i^) then there is a 

constant eo for each < e < eo, the problem i9.6]) has a unique solution /2 € W'^'''{T) with the 
property 

ll/2||L~(r) + II /2||L-(r) + II /2|U,(r) < C\\h\\L,^r)- 

Proof. Under the non-degeneracy condition of F in (|1.3p , the existence and the a priori estimates 
can be easily proved. □ 

We then consider the following problem 

Be 1 

£*^^{e) = -e^A^e - Xqc = g in F, — + -le = on dT. (9.7) 

OT 2 
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Lemma 9.2. If g (z i^(r) then there exists Eq > such that for a sequence {si)i with < £; < Eq 
such that problem fff. 7| ) has a unique solution e € W'^''^{T) which satisfies 

\\e\\b < C e-^ ||ff||L<,(r)- 

Moreover, if g e W^''^{T) then 

e?ll^^e||i,,(r) + ei||V^e||i,(r) + ||e||Loo(r) < C\\g\\w2,^r)- (9.8) 
The proof will be given in subsection 19.21 □ 

We finally consider the following system 



C*{f2,e)^{Cl{f), CUe)) ^{h,g) in T 
^+//2-S, |£ + i/e = onar, 

OT OT 2 



(9.9) 



where ^ is a smooth function. 



Lemma 9.3. Under the nondegeneracy condition il.3]} . if h, g € i'^(r) then for the sequence of 
the parameter e in Lemma \9.2[ there is a unique solution (/2, e) in W'^'''{T) to problem i9.9\) which 
satisfies 

II/2IU + llellfc < C[ \\h\\L,ir)+s-'\\g\\LHr) + llS|U,(r) ]. 
Proof. Under the non-degeneracy condition (|1.3p . there exist /2 and eg satisfying 

A^/2 = 0, 1^+7/2 = S. 

Setting /2 = /2 + /o to the system (|9.9p . the final conclusion can be derived from Lemma I^TTl and 
Lemma EH □ 



Proof of Theorem ll.lt Let (/, e) S F, where F is defined in p.38p . and define 

( M/2,e), 5(/2,e) ) ^ [sAuif2,e)+eKu{f2,i), e^A2,{f2,e)+e^K2eif2,S) ), 
S = $n(/2,e~). 

From (j9.5p . Ai^ and are contraction mappings of its arguments in F ■ By Banach Contraction 
Mapping theorem and Lemma 19.31 we can solve the nonlinear problem 

C*if,e)^(^CUf2),C*2M)^ih,g), 

with the boundary conditions defined in (|9.9p on the region F- Hence, we can define a new operator 
Z from F into F by Z(/2,e) — (/2,e). Finding a solution to the problem (j9.ip - (|9.4p is equivalent 
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to locating a fixed point of Z. Schauder's fixed point theorem applies to finish the proof of its 
existence. Hence, by Proposition 17. II and the arguments followed, we complete the existence part 
of Theorem ll.il Other properties of in Theorem 1 1 . 1 1 can be showed easily. □ 



9.2 Proof of Lemma lOl 



To prove Lemma 19.21 we follow the method introduced in |12j , which relies only on elementary 
considerations on the variational characterization of the eigenvalues of the operator /Zje ^-i^d the 
Weyl's asymptotic formula in (|2.20p . We remark this approach is a slightly different from [24 and 
[25] where Kato's theorems were the main tools. 

First, we consider the following eigenvalue problem 



dv 1 

Xqv = A^v in r, + -/w 

OT 2 



on dr. 



(9.10) 



We denote its eigenvalues Ag j- in non-decreasing order and counting them with multiplicity. Here 
Aq is the unique positive eigenvalue to the eigenvalue problem (jl.Sp . which implies the spectrum of 
contains negative or zero eigenvalues. From the Courant-Fisher characterization we can write 
Agj in two different ways: 



A, 



sup 



inf 



2e' 



(9.11) 



Agj- — _inf 

Here 2j (resp. S^-i) represents the family of j dimensional (resp. j — 1 dimensional) subspaces of 
_ff^(r) constituting of functions defined on F with boundary condition in (j9.10p . and the symbol 
_L denotes orthogonality with respect to the scalar product. There holds the following result 
for the estimates of gap between two successive eigenvalues. 

Lemma 9.4. There exits a number Sq > such that for all < ei < £2 < o,nd all j > 1 the 
following estimate holds. 




(9.13) 



In particular, the functions e G (0, Eg) ^ ^e,j 0,'^^ continuous and increasing. 

Proof. Let us consider small numbers < ei < £2- We observe that for any v with Jp = 1, we 
have 

^2^ ^^*2e,y ~ y^2e,V = Ao(£r' - ^2')- 



sup 

t) G S, v=iO 



It ^^2e^ 



(9.12) 
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Then the resuh follows. □ 

Proof of Lemma I9.2t For £ £ N, choose ai = 2~^. In order to find a sequence of values 
G ((T£+i,CT^) such that the spectrums of the operators C2^^, for large i, stay away from 0, we 
define 

JC] ^ {e £ {ai+i,ai) : ker/:2e 7^ ^ }: {<7e+i,cri) \ JC\. 

It is crucial to estimate the cardinality oi JC\. If e G IC] then for some j we have that = 0. 
The monotonicity of the function e G (0,eo) ^ A^j implies that Ao-f^u < 0. Hence, 

cwA{lC\) < N,,^,, (9.14) 

where is the number of negative eigenvalues of the operator /Ije- 

We now give an asymptotic estimate on the number of negative eigenvalues of the differential 
operator C^^. By {pi)i we will denote the set of eigenvalues of the eigenvalue problem 

-A^cj = pLu in r, + -/w on dV. 

OT 2 

From the Weyl asymptotic formula as those in (|2.20p and the formula in l|9.12p . one derives 

N, > Cr(l + o(l))e-2, 

where Cr is a fixed constant depending on the volume of the manifold F and its dimension. To 
prove a similar upper bound, we choose i to be the first index such that e'^pi — Aq > 0. Then from 
the Weyl formula we find that 

i = Cr{l + o{l))e~^. 

Define = spanjwf : Z = 1, 2, • • • , j — 1}. For an arbitrary function v G iJ^(F) and v _L Sj-i, 
we can write 

W = ^ KiUJi. 
l>3 

Plugging this v into (|9.1ip and using the Weyl formula, we also have 

Ne < Cr(l + o(l))£-^ 

Hence we get that 

~ Cre"^ as £ ^ 0. 

The last inequality and (|9.14p imply that card(/C^) < Ccr^^, and hence there exists an interval 
(flf, bi) such that 

{a,M)clCj, \b,-ai\>'^^^^^>2C^al (9.15) 

card(AL£) 



49 



for a universal positive constant Co, independent of t. By setting si — {ai + bi)/2 for all large 
£ e N, we conclude that C^a invertible and there exists a number C > 0, independent of £, such 
that for all j E N there holds 

|Ae,j| > Cej. (9.16) 
Assume the opposite, namely that for some j we have 

|Ae,,j| < Sej, 

with 5 arbitrarily small. Since S K-f, then |Ae^.j| > 0. Let us assume that 

< Ae,,j < 6ej. (9.17) 

Then from Lemma |9.4[ we have 

Aa„, = A,,,, - M-^(A,„, + Ao). (9.18) 
The inequalities in (j9.15p and (|9.17l) imply that 



A.,,, < fc| - Coa|^^^^^(A,,, + Ao) < 0, 

if 6 is chosen a priori sufficiently small. It follows from the continuity of the function e e (0, eo) 
Agj that must vanish at some e G (a£,££), and we get a contradiction with the choice of the 
interval (ae, bg). 
The case 

-Saj < Kei,] < 
can be handled similarly. In fact, we have the inequality 

K,j = A,,j + ~/'^ (Ae„j- + Ao) > 0. 

Hence, the proof of (|9.16p for the spectral gap between critical eigenvalues was complete. 
As a consequence, the solution to (|9.2I) exists and satisfies 

||e||L.(r) < Cei^MLHD- (9.19) 
From (I9.19|) by a standard elliptic argument one can show 

Si^ \\A^e\\mr) + £i\\y^e\\LHr) + ML^iT) < C^e^" ll^alU^cr)- 
The reader can refer [12] for proof of further estimate in (|9.8p . □ 
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A A Linear Model Problem I 



Recall that w is the even function defined in (jl.4p and Z is the even eigenfunction defined in the 
eigenvalue problem (|1.5p . Recall that, A represents the strip in M"^ with the notations 



A = {{x,0,r)) : X eR, {0, ry) G 9ID) x M+ }, 
dA ^ { {x, 9,7]) : X eM., 6 e dB, f] ^0 }. 

We first consider the following linear problem 



(A.1) 



PL + j^^ee + - K^" + Pii'^- V° =0 in A, (A.2) 

(/>° = G{x, 9) on dA, (A.3) 
(f>°{x,9,T])w^{x)dx ^0, [ (j)°{x,9,r])Z{x)dx ^0 in T^, (A.4) 



where > Ao + 1 is a large positive constant, li{d) and hiO) are two smooth positive functions 
given in (I4.3p . Suppose the following orthogonality conditions hold 

G{x,0)w^{x)dx = 0, J G{x,0)Z{x)dx ^0. (A.5) 

Lemma A. 1.1. If G E L^{dA), and the orthogonality conditions iA.5\) hold, then there is a 
unique solution (jp to the problem iA.^] - ^.^^ for any large positive constant K . Moreover there 
is a constant C > 0, independent of e, such that the solution to the problem \A.2\l - lA^ satisfies 
a priori estimate 

ll/llff^(A) <C||G|U.(aA). 

Proof. Since K is large, the proof of the existence and uniqueness of the solution to (jA.2p - (IA.4p 
and its estimate is standard. To show the L^-orthogonality (|A.4p . using the equations of Z[x) and 
(jp and also the condition ()A.5[) . for 

^(^,77)= / 4P{x,9,7])Z{x)dx, 

one finds 

-j^Vee + -j^^'frir, - {K - 1 - Xo) ip = inT,, ip^{9,0)^0. 
Choosing X > Aq + 1, we deduce that 



ip{9,r]) = / (l)"{x,e,T])Z{x)dx = in 



51 



Similarly we have 

(j)^{x,9,ri)wx{x)dx = in Fg. 



□ 

A special case of Lemma I A. 1 . 1 1 is the following problem: finding function £ _ff^(A) such that 

1-1 

-(pri-q ^ i^<P + pw^ ^0=0 in A, 

(A.6) 



(j)^ = G{x, 6) on 9A, 

where if is a large positive constant. 

Lemma A. 1.2. Suppose the function G{x,9) is even in the variable x, then there exists a large 
positive constant K such that the problem \A.9j) has a unique solution (j), which is an even function 
in the variable x and satisfies 

||0||//2(A) < C||G||l2(9a), 
Moreover, if G{x, 0) is exponentially decaying in x, then 

|0(a;,6',7?)| < Ce-"l^l, (A.7) 

where a > and the constant C does not depend on e. 

Proof. For any K large enough, the proof of the existence and uniqueness of the solution to (jA.6| 
and its estimate is standard. By uniqueness and evenness of G{x, 6), is an even function in the 
variable x. By the exponentially decaying of G{x,9), we also have (|A.7p . □ 
Next, we consider the following problem 

Co{4>) = (j)xx + 727HV'^««+727^'^w ~ + Pw^~^(i) = h in A, 

4>^ = G{x, 9) on aA, (A.8) 
'{x,9,vi)wx{x)(ix — I (l){x,9,ri)Z{x)dx — in F^. 



Lemma A. 1.3. If h E L^(A), G E L^{dA) and the orthogonality conditions HA . 5]) hold, then for 
any solution (j) to problem \A.^) we have 

II^IIh^(A)<C[ ||/l||L2(A) + l|G'|U2(aA) ] 

where the constant C does not depend on h, G and e. Furthermore, if \h\ + \G\ < Ce~"'^', then 
\4>\ < Ce-="l^l for some C, c> 0. 
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Proof. Let (jp(x,9,ini) be defined in Lemma lA. 1.11 and (p = (fp + (f>. Then we have 

'1'^^ + 727ST'?^ee + ■m^'Pvv ^ ^ + pw^^^4> = ^ + (1 - K)(ff in A, 

(j)rj =0 on dA, (A.9) 
J (f>{x,9,ri)wx{x)dx — 0, J (f>{x,9,ri)Z{x)dx = in Fg. 

Let £_k,k — 1,2,--- , be the eigenfunctions (corresponding eigenvalues Xk,k = 1,2,---) of the 
following eigenvalue problem 

Let us consider Fourier series decompositions for h and (p of the form Let us consider Fourier series 
decompositions for h + {1 — K)(jP and <p: 

oo 

(l){x,0,r]) = ^(I)k{x)(,k{e0,er]), 

oo 

h{x, e,Tj) + {l- K)4P{x, e,ri)^Y. hk{x)(,k{ee, STj). 
From the equation (jA.9[) we arrive at the following equation 

- s^Xk(j)k + 4>k,xx - 4'k+ pw^^^(i)k = hk, (A.IO) 
with the orthogonality condition 

(f>k{x)wx{x)dx — 0, / 0fc(x)Z(x)da; = 0. (A. 11) 

: Js. 

Let us consider the bilinear form in iJ^(R) 

Since (|A.11[) holds uniformly in k we conclude that 

C[ \\4>k\\l2(^K) + ||<?!>fe,x||i2(R) ] < B{(j)k,(j)k), 

for a constant C > independent of k. Using this fact and equation (jA.lOp we arrive at 

(1 + Xle^)Uk\\hm + UkAhiw) < C\\hk\\h(^y (A.12) 
Moreover, we see from (jA.10|) that (pk satisfies an equation of the form 
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where ||/i/c||l2(r) < C||'ifc||L2(R)- Hence it follows that 

II 2:3: II i2(R) < C||^fc|li2(R)- (A. 13) 

Summing up estimates (IA.12p and (|A.13I) in k, we conclude that 

P'0lli^(A) + ll^'/'lli^(A) + ll'/'lli^(A) < cmi^^^y 

The final estimate follows from the estimates of (p and 0" . □ 
A corollary of Lemma IA.1.31 is the following 

Corollary A. 1.1. Let G G L^{dA) satisfy the orthogonality conditions HA . 5\) and h — 0. Then 
problem SA.8\) has a unique solution 4> such that 

|l0llff2(A) < C||G'||i2(9A), 

where the constant C does not depend on G and e. Furthermore, if \G\ < Ce~"'^', then \(f)\ < 
(-g-cal^l gg^g C, 0. □ 

Appendix B. A Linear Model Problem II 

Recall & represents the strip 

{ {x,z) : X eR,z e Te} (B.l) 

in M^. 96 is the component of the boundary of (3, i.e. 

de ^ { {x,z) : X eR, z e dVe}. 

We consider the following problem: given h e L'?(6) and G £ L'^{d&), finding functions 
(j) G W^^i{e), c,de Li{r) and Ai, A2 such that 

+ (f) pw'P^^ (j) ^ h + c{eri) xi^l^l) + d{ef]) x(e|a;|) Z in 6, 

1^ = G on de, (B.2) 

/ (t>{x, z)wxdx — Ai, / (j){x, z)Z{x)dx = A2 in Fe, 
Jr Jr 

where denotes the inward normal of d&, and x(t) a smooth cut-off function such that xit) = 1 
for |t| < 10a and %(<) = for t > 20a, and cr > is a small constant defined in Section [5l 

Lemma B.2. 4. There exist functions c{ez), d{sz) with respect to h such that the problem HB. 2^ 
has a unique solution (j) = Ti{h, G). Moreover, 

MWe < C^(l|/i||L.(e) + ll5llL.(a6)), 
||Ai||L9(r,) < C'(ll^llL'!(e) + ll5llL'j(a6)), Vj = l,2, 
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where the constant C does not depend on h, G and e. 
Proof. The proof will be carried out in three steps. 

Step 1: Let us assume that in problem (jB.2p the terms G, c{£z),d{ez) are identically zero. Arguing 
as in Lemma [A. 1.31 for any h e L'^{&) and any solution </) e W'^''^{&) of problem (|B.2[) we have 

MW,<C\\h\\L.ie) (B.3) 

Step 2: We claim that the a priori estimate obtained in Step 1 is in reality valid for the full 
problem (IB.2|) . We first choose suitable Ai and A2 such that 

V^^Ai = / g{x,z)w^dx, V^^Aa = / g{x,z)Z{x)Ax in T^. (B.4) 
■/r Jr 

Let (fP be the solution of 



Note that we have 



6°||,,,<C||G|U,(ae) (B.5) 



By defining 

= / (fP [x^ z)wx<Ax^ ^2— (t)^ {x , z) Z {x)dx in Fe, 
Jr Jr 

we have that, to prove the general case it suffices to apply the argument with 

- (Ai - Ai)wx (A2 - A2)Z{x) 

^ = ^-^ + /,.^d. + J,zHx)dx ■ 

The (p satisfies a problem of a similar form with homogeneous Neumann boundary condition and 
orthogonality condition, as well as h replaced by a function h with norm bounded by 



\Li{e) < C* (||ft-||L9(6) + \ \G\\LHde))- 
Step 3: We consider the problem 

(f>xx + A'"^^ — (j) + pw^~^4' = h + cxWx + dxZ in ©, 
— — =0 on o&, 

{x, z)wx<ix — 0, / (t>{x, z)Z{x)dx ^ in Fg. 
Jr 

The existence of the solutions can be proved similarly as that in Lemma [A. 1.31 There also hold 
the priori estimate 

mks<C{\\h\\L.(e) + \\9\\LHde)), 
||Ai||L<,(r,) < ^(11/111^,(6) + WgWLUde)), = 1,2. 

□ 
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